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^ ' Abstract 



A detailed error analysis is carried out for the determination of the MS charm 
quark mass mc(mc) from moments at order of the charm cross section in 
e'^e~ annihilation. To estimate the theoretical uncertainties the renormaliza- 
tion scale is implemented in various ways including energy-dependent func- 
tions, which lead to "contour-improved" predictions. We obtain rrlcijnc) = 
1.29 lb 0.07 GeV which contains a substantial theoretical uncertainty. 



Introduction 



At the ongoing and future B-physics experiments a realistic estimate of the uncer- 
tainties in the values of the bottom and charm quark masses will become increasingly 
important for the measurements of the CKM parameters and the search for new 
physics, particularly from inclusive B-decay rates. ^1 12] Moments of the respective 
e+e" R-ratios 0111, 



where Rqq = cr(e"'"e~ qq + X)/cr(e"*"e~ are the most important instru- 

ment at present to determine the bottom and charm quark mass parameters. (See 
also Refs. IH] for recent bottom quark mass extractions from spectral moments in 
semileptonic B-decays.) 

In general, one can distinguish between two regions in n, which require a different 
theoretical treatment. For low values of n the moments are dominated by relativistic 
dynamics and scales of order of the heavy quark mass rriq. This allows the use 
of the usual expansion in the number of loops for the theoretical computations, 
and the MS scheme is an appropriate choice for the heavy quark mass parameter. 
However, the lack of data for Rqq in the continuum regions above the quarkonium 
resonances introduces model-dependent errors that have not been quantified in most 
of the previous analyses. (See Refs. for recent reviews.) On the other hand, for 
large values of n the continuum regions are suppressed and the moments become 
dominated by the quarkonium resonance region where good sets of data have been 
obtained in the past. However, the theoretical predictions of moments for large 
values of n is more complicated since the usual loop expansion breaks down and 
the size of non-perturbative effects increases. Here, summations of higher order 
contributions proportional to powers of (a^-y/n) need to be carried out in order to 
capture the relevant non-relativistic perturbative information |Hl E] , and so-called 
threshold masses fOl ^1 ^] are appropriate schemes for the heavy quark mass 
parameter. Moreover, there is an upper "duality" bound for the possible choices 
of n since the energy range contributing to the moments, which is of order mq/n, 
needs to be larger than the typical hadronization scale Aqcd- P UH] In the case 
of Rbb and the determination of the bottom mass, this bound is around n = 10, 
and the two ranges in n are believed to be well separated with their boundary 
being approximately at n = 4. A good number of analyses exists for small and large 
values of n and respecting the cancellation of the dominant renormalon contributions 
associated to the choice of the quark mass definition. [TH [T3 CHI HZl d] (See also 
Ref. [in] for an early analysis using the MS bottom quark mass for larger values of 



For Rcc and the determination of the charm mass, the distinction between large- 
and low-n moments is more delicate because rric is not much larger than Aqcd- 
Here, the upper duality bound for n is around 3 or 4. This leaves basically no 
space at all to carry out the non-relativistic summations that can be applied in the 
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bottom quark case because the corresponding techniques are only vahd for large n 
and as long as rug/n is larger than Aqcd- On the other hand, even for n < 4 the 
non-relativistic region close to the cc threshold can have a considerable contribution 
to the moments, while the model-dependences from the experimentally unknown 
continuum region can still be significant. For a reliable (error) analysis these issues 
need to be taken into account. 

In this paper we determine the MS charm quark mass fnd^c) from moments 
Pn with n < 4 using the usual loop expansion in powers of as- In a previous 
work an error of less than 30 MeV was obtained using results at order 
and a fixed renormalization scale /i between rric and Arric. This analysis assumed 
that the moments are dominated by scales of order m^. In our analysis we focus 
on the theoretical uncertainties coming from the low-energy region close to the cc 
threshold. It is well known that close to threshold, where the cm. velocity (3 of 
the quarks is small, the dominant perturbative contributions to i?cc are in fact 
governed by renormalization scales of the order of the relative quark momentum, 
/i ~ rricP. [201 1^ We believe that this effect should not be ignored in determinations 
of the charm quark mass since it accounts for important higher-order information. It 
is the main aim of this work to study the impact of these higher-order contributions. 
In particular, we find that they affect the estimate of the theoretical error. We also 
carry out a careful study of all sources of uncertainties including the continuum 
regions where no data is available. 



Theoretical Moments 

For the QCD parameters used in this work we adopt the MS renormalization scheme 
and the convention that the charm quark participates in the running (n/ = 4). The 
theoretical moments are directly parametrized in terms of Trie = rnc{mc), and the 
masses of the up, down and strange quarks are set to zero. For the running of the 
strong coupling we use four-loop renormalization group equations, and three-loop 
matching conditions at the scale mfe(mf,) = 4.2 GeV where we switch to rij = 4 
flavors. We use three different methods to implement the renormalization scale 
for the theoretical predictions of the moments and apply order results obtained 
earlier in Refs. [201 1221 123 121 • 

The first method is simply based on a fixed energy-independent renormalization 
scale, /i^ ~ ml, in analogy to previous analyses. In the OPE including the first non- 
perturbative contribution from the dimension-four gluon condensate, the moments 
have the form 
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n 


1 


2 


3 


4 


J n 


1.0667 


0.4571 


0.2709 


0.1847 


J n 


2.5547 


1.1096 


0.5194 


0.2031 


J n 


2.5896 


2.7790 


1.6390 


0.7956 


9" 


-16.042 


-26.737 


-38.890 


-52.352 


9i 


-15.028 


13.006 


78.710 


190.935 



Table 1: Coefficients of the theoretical expressions for the moments Pn at order a1 
for massless light quarks based on a fixed energy-independent renormalization scale. 
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Here, (3o = 11 — 2/3n/. The order as terms were given in Refs. [SI E] and the 
coefficients for the condensate contribution were taken from Refs. |25[ OHj adopting 
the renormalization-group-invariant normalization of the gluon condensate. The 
numerical results for the coefficients in Eq. Q forn = 1, 2, 3, 4 are collected in Tab.^ 
The numbers for /° and agree with results given in Ref. ^1]. The numbers for 
differ because we also include the contributions from secondary cc production, where 
the charm pair is produced through gluon splitting off primary massless quarks. [221 
For the determination of mc(mc) described in subsequent sections we use /x^ = ^^M^ 
with M = 1.3 GeV as the renormalization scale and vary the parameter ^ in order 
to estimate the perturbative uncertainties. 

The second method is based on an energy-dependent renormalization scale of 
order of the cm. energy, /i^ ~ s. We note that the same energy-dependent scale 
also needs to be used for the implementation of the MS charm mass mcirfic) in order 
to avoid an incomplete cancellation of the large higher order corrections associated to 
the (!?(Aqcd) renormalon of the pole mass definition. A technical issue is related to 
the fact that in the MS scheme for the charm mass parameter, Rcc is more singular at 
the threshold s = Arfic than in the pole scheme. This makes the computation of the 
perturbative contributions of the moments based on the dispersion integral in Eq. (|T|) 
somewhat cumbersome. It is therefore advantageous to carry out the computation 
using a contour integration in the positive complex s-plane (Re(s) > 0) around the 
cut of the vacuum polarization function IIcc associated to Rcc, 

„ 1 r ds 



1 n,,, _ 

2Tci Jc V4m 



^ ^ ^ n,,(.-). (5) 



(4m2)n 2TTi Jc s' 

The same procedure can of course also be applied for the ffist method. The path C 
of the integration is illustrated in Fig.Q For practical purposes we use s = s/{4mc) 
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Figure 1: Path of integration in the complex s-plane for the computation of the 
moments. 



n 


1 


2 


3 


4 




0.2433 (0.1839) 


0.1093 (0.0822) 


0.0498 (0.0379) 


0.0167 (0.0135) 


CI' 


0.0434 (0.0223) 


0.0308 (0.0169) 


0.0165 (0.0094) 


0.0071 (0.0042) 




0.0235 (0.0134) 


0.0107 (0.0061) 


0.0047 (0.0027) 


0.0012 (0.0009) 




0.0198 (0.0116) 


0.0082 (0.0046) 


0.0044 (0.0024) 


0.0027 (0.0015) 




0.3105 (0.2183) 


0.2339 (0.1400) 


0.1917 (0.0968) 


0.1654 (0.0685) 




0.0385 (0.0190) 


0.0658 (0.0283) 


0.0759 (0.0270) 


0.0820 (0.0239) 




0.0377 (0.0187) 


0.0405 (0.0158) 


0.0428 (0.0136) 


0.0455 (0.0117) 




0.0000 (0.0000) 


-0.0377 (-0.0187) 


-0.0593 (-0.0252) 


-0.0756 (-0.0277) 



Table 2: Coefficients of the theoretical expressions for the moments P„ at order 
ct^ for massless light quarks based on energy-dependent renormalization scales as 
described in the text for ^ = 1 (^ = 2). 



as the integration variable and ^ = ^'^M'^s as the renormalization scale. This 
choice allows to explicitly factor out the dependence of the moments on rfic- The 
perturbative part of the moments can then be written in the form 

.^2^2 
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and 0{a1), respectively. Exem- 



(i.e. a7 (4.2GeV) = 0.225), the 



where the superscripts 0, 1,2 refer to 0(1), 0{as 
plarily, for ^ = 1 (^ = 2) and a7^^{M,) = 0.118 
coefficients are given numerically in Tab.|21 

The third method is based on an energy-dependent renormalization scale of 
order of the relative three- momentum of the charm quarks, fi^ ~ —s/S"^, where 
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(3 = ^1 — Ami/ s. Again, the perturbative part of the moments is computed using 
the contour integration in Eq. (jH}. In the complex /i^-plane, ^^(/i^) has a cut along 
the negative real axis. It is therefore necessary to introduce the negative sign. With 
this choice the cut in as(yU^ ~ — agrees with the one of the vacuum polarization 
function, and contours can be found such that (3 is never small along the path, 
and the perturbative description is appropriate.^ In fact, this is in close analogy to 
the so-called "contour-improved" approach to compute the hadronic r-decay width 
(see e.g. Ref. |2Z1 EHl)- For our numerical computations we use s = s/(4m^) as 
the integration variable and fi^ = ^^M^(l — s) as the renormalization scale. The 
perturbative part of the moments can then be written in the form 

;^f + ^(c^fln(^)+Cf)]}. (7) 

Exemplarily, for ^ = 1 = 2) and a"-^=^(M,) = 0.118 (i.e. 0^-^=^(4.2 GeV) = 
0.225), the coefficients are given in Tab. El Note that CP = for any choice of the 
parameters because the corresponding integrand does not have a pole at s = due 
to the additional factor of ln(l — s). 
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Experimental Moments 

For the contributions to the experimental moments from the J/ if) and the we use 
the most recent averages for masses and e^e~ widths given by the PDG {Mj/^ = 
3.0969 GeV, M^, = 3.6860 GeV, 1})^" = 5.14 ±0.31 keV, T^t^" = 2.12 ±0.12 keV). 

For the energy region between 3.73 and 4.8 GeV we employ the data for the total 
R ratio, -Rtot? obtained by the BES collaboration. 30 To obtain experimental data 
for Rcc we use the last four data points below 3.73 GeV to fit the non-charm R ratio, 
Rnc, assuming it to be energy- independent, and subtract it from the numbers for 
Rtot- The effects from the Z exchange and the small logarithmic energy- dependence 
expected from theory are below 1% and can be neglected. The determination of the 
statistical uncertainties of Rcc is standard. The determination of the systematical 
uncertainties in Rcc is more subtle because only systematical errors for _Rtot (with 
all individual contributions being added quadratically) but no specific numbers for 
the D-meson final states above 3.7 GeV were given in Ref. jHOj. Thus, naively sub- 
tracting the systematical errors of i?nc from the ones in i?tot underestimates the 
actual systematical errors in R^c- To obtain numbers for the systematical errors in 
Rcc we make the assumption that the systematical errors for the non-charm and 
the charm final states are equally large, i.e. the systematical errors from Rcc are 
obtained by multiplying those of -Rtot by a factor 1/a/2. This procedure is supported 
by the observation that the relative systematical errors given in Ref. [30] for energies 

^For the first and the third method one can also choose a closed path around the origin in the 
counter-clockwise direction. 
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contribution 


Pi 

X 102 GeV^ 


P2 

X 10^ GeV^ 


X 10^ GeV^ 


P4 

X 10^ GeV^ 




8.82(38|38) 


9.20(39|39) 


9.59(41|41) 


10.00(43|43) 


^' 


2.16(9|9) 


1.59(6|6) 


1.17(5|5) 


0.86(3|3) 


3.7-4.8 GeV 


3.23(6|33) 


1.80(3|19) 


1.02(2|11) 


0.59(1|7) 


4.8 - 7.25 GeV 


3.82(25) 


1.19(8) 


0.39(3) 


0.13(1) 


7.25 - 10.5 GeV 


1.45(3|10) 


0.21(0|1) 


0.03(0|0) 


0.00(0|0) 


10.5 GeV - Mz 


1.27(1) 


0.06(0) 


0.00(0) 


0.00(0) 


Mz -oo 


0.02(0) 


0.00(0) 


0.00(0) 


0.00(0) 



Table 3: Individual contributions to the experimental moments including statisti- 
cal and systematical uncertainties. If only one error is given, it comes from theory. 
The resonance contributions have been determined in the narrow width approxima- 
tion, {Pn)k = QTrr^^*^ /[a^(3.1 GeV) M^""*""^], where for the electromagnetic coupling 
[a(3.1 GeV)]-i = 134.3 has been adopted. 



above 3.7 GeV are in average approximately a factor 1.5 larger than expected from 
extrapolating the relative systematical errors for energies below 3.7 GeV. 

In the region above 4.8 GeV the only other data available for i?tot is from the 
MDl experiment [HI] and from the CLEO collaboration [32] in the energy region 
between 7.25 and 10.5 GeV. For the determination of i?cc and the uncertainties we 
use the same method as described in the previous paragraph. The resulting relative 
systematical errors amount to about 8%. 

There is no continuous experimental data for i?cc in the regions between 4.8 and 
7.25 GeV and above 10.5 GeV. As shown below, these contribution turn out to be 
not very important. For our analysis we use a model for these contributions based 
on the theory prediction of Rcc at order a^. This is reasonable since it is known 
that perturbative QCD agrees with data for cc production to several percent at Mz 
and at the level of a few 10% percent in the LEP2 region. 33 . The resulting model- 
dependence in the determination of is negligible for n > 1, as shown below. 
For n = 1 the contribution from the region between 4.8 and 7.25 GeV leads to a 
model-dependent error of about 10 MeV. 

In Tab. IS] a collection of all contributions to the first four moments is given. For 
contributions coming from experimental data the first error is statistical and the 
second systematical. For the contributions from the J/tp and the ip' we used the 
averaged PDG [21] errors. Since detailed numbers for the various error sources are 
difficult to obtain, we make the simplified assumption that statistical errors (be- 
ing later combined in quadrature) and systematical errors (being later combined 
linearly) are equal. This is consistent with information given in the original litera- 
ture. For the continuum regions that have to be estimated from our theory-model 
the contributions from below and above Mz are distinguished in order to visualize 
the impact of these two regions. The respective errors given in Tab. [HI come from 
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n 


1 


2 


3 


4 


Method 1 (3/4 < ^ < 3) 


1.286(100) 


1.280(66) 


1.274(51) 


1.269(44) 


Method 2 (3/4 < ^ < 3) 


1.280(101) 


1.272(61) 


1.267(47) 


1.264(41) 


Method 3 (3/4 < ^ < 3) 


1.287(100) 


1.296(69) 


1.298(58) 


1.298(53) 


combined (3/4 < ^ < 3) 


1.283(104) 


1.288(77) 


1.288(68) 


1.287(64) 


Method 1 (1 < C < 4) 


1.277(94) 


1.272(59) 


1.267(47) 


1.264(45) 


Method 2 (1 < ^ < 4) 


1.268(92) 


1.264(56) 


1.262(45) 


1.260(43) 


Method 3 (1 < ^ < 4) 


1.278(94) 


1.291(65) 


1.294(55) 


1.295(52) 


combined (1 < ^ < 4) 


1.274(98) 


1.282(74) 


1.283(66) 


1.282(65) 



Table 4: Results for md^c) for different methods to implement the renormalization 
scale. The uncertainties represent the sum of all individual error sources. The 
combined result represents the total mass range covered by the three methods. All 
numbers are given in units of GeV. 



varying the strong coupling in the range as{Mz) = 0.118 ± 0.003, the MS charm 
mass in the conservative range nic = 1.3 ± 0.2 GeV and the renormalization scale 
between /x^ = s/16 and /i^ = s. The central values are obtained from the average 
of the corresponding extremal values. The relative error obtained in this way is not 
larger than 7% for all three regions where theory is employed. Note that for our 
error estimates for Trie (see Tab.E} we adopt, by hand, a relative error of 10%. This 
reduces the mo del- dependence of the n = 1 results to an acceptable level. 

Uncertainties in mc{rric) 

For the determination of fnc{fnc) we fit single moments using the three different 
methods for the theoretical predictions described above. The contribution from the 
gluon condensate is taken from the first method for all cases. We use as^~ [M^) = 
0.118±0.003 (as^=^(4.2 GeV) = 0.225±0.012) and (^G^) = (0.024±0.024) GeV^ as 
theoretical input. Our final results for n = 1, 2, 3, 4, for the three different theoretical 
methods and for varying the renormalization scale parameter in the ranges 3/4 < 
,^ < 4 and 1 < ^ < 4: are given in Tab.|31 The central values are obtained using 
the respective central values for all experimental data and the theoretical input 
parameters, and by taking the median of the range in Wic^jn^ covered by varying ^ 
in the respective ranges. The uncertainties are obtained by combining all individual 
error sources as described below and shown in Tab. El In Tab.|3]we have also given 
combined results representing the total range in rfic covered by the three theoretical 
methods. 

A detailed account of individual errors is presented in Tab. El For experimen- 
tal data the first error is statistical and the second systematical. The results are 
only given for the second and the third theoretical method, as the first and second 
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Method 2 


Method 3 


n 


1 


2 


3 


4 


1 


2 


3 


4 


3/4 < ^ < 3 


central 


1280 


1272 


1267 


1264 


1287 


1296 


1298 


1298 




12/12 


9/9 


7/7 


6/6 


11/11 


9/9 


7/7 


6/6 


ip' 


3/3 


1/1 


1/1 


0/0 


3/3 


1/1 


1/1 


0/0 


3.73 - 4.8 GeV 


2/10 


1/4 


0/2 


0/1 


2/10 


1/4 


0/2 


0/1 


4.8 - 7.25 GeV 


12 


3 


1 





11 


3 


1 





7.25 - 10.5 GeV 


1/3 


0/0 


0/0 


0/0 


1/3 


0/0 


0/0 


0/0 


10.5 GeV - oo 


4 











4 













4 


8 


11 


14 


4 


6 


8 


9 


Sas{Mz) 


14 


8 


5 


3 


17 


17 


17 


18 




28 


18 


12 


9 


25 


19 


15 


12 


combined error 


101 


61 


47 


41 


100 


69 


58 


53 


1 < ^ < 4 


central 


1268 


1264 


1262 


1260 


1278 


1291 


1294 


1295 




12/12 


9/9 


7/7 


6/6 


12/12 


9/9 


7/7 


6/6 




3/3 


1/1 


1/1 


0/0 


3/3 


1/1 


1/1 


0/0 


3.73 - 4.8 GeV 


2/10 


1/4 


0/2 


0/1 


2/10 


1/4 


0/2 


0/1 


4.8 - 7.25 GeV 


12 


3 


1 





12 


3 


1 





7.25 - 10.5 GeV 


1/3 


0/0 


0/0 


0/0 


1/3 


0/0 


0/0 


0/0 


10.5 GeV -oo 


4 











4 













5 


9 


14 


20 


5 


7 


10 


14 


Sas{Mz) 


11 


5 


3 


2 


11 


10 


9 


9 




24 


15 


10 


7 


24 


20 


18 


14 


combined error 


92 


56 


45 


43 


94 


65 


55 


52 



Table 5: Central values and individual errors for mc(mc) in units of MeV based on 
the methods described in the text. 



methods yield very similar results, particularly for n > 1. This is expected since 
the high-energy range, s 3> 4m^, where the renormalization scales of the two meth- 
ods differ, is suppressed. The uncertainties from experimental data and from the 
treatment of the non-perturbative corrections are similar for the three theoretical 
methods. As mentioned before, the model-dependent error coming from the lack 
of experimental data in the continuum region is negligible for n > 1 and amounts 
to only about 10 MeV for n = 1. This is in contrast to determinations of Tnb{rni,) 
where the model-dependent error is significantly larger, particularly for n = 1 and 
n = 2. |T3] The uncertainties coming from the error in as{Mz) for n > 1 are larger 
for the third method than for the first and the second one, since variations in the 
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value of as{Mz) have a larger impact at lower scales. To obtain combined errors we 
have added all statistical uncertainties in quadrature, and all other uncertainties lin- 
early. Since the statistical errors are smaller than the other sources of uncertainties, 
the reader should note that our total errors cannot be interpreted statistically. 

It is instructive to have a closer look onto the perturbative uncertainties obtained 
by the different theoretical methods. For all methods the uncertainties obtained by 
varying ^ either in the range 3/4<^ <3orl<^ <4 become smaller for larger 
values of n. This behavior is, however, less pronounced for the third method. Thus, 
each of the methods yields a behavior in contrast to the expectation that the scaling 
uncertainties grow like the uncertainties from the non-perturbative corrections, since 
for larger n the moments become dominated by dynamics at lower energies. On 
the other hand, comparing the results for the central values obtained by the three 
theoretical methods it is found that they tend to decrease for larger values of n 
for the first and second method, while there is the opposite effect for the third 
method. In fact, the difference between the central values for n > 1 is larger than the 
uncertainties based on varying ^. This feature is related to the different treatment 
of higher-order contributions and appears to be similar to the discrepancies found 
between the "contour-improved" and the "fixed-scale" approaches for predictions 
of the hadronic r-decay rate. (See e.g. Refs. jSH EH] for discussions on this issue). 
Interestingly, the convergence behavior of the moments for all three methods is quite 
good. For example, for n = 3, nic = 1-3 GeV, ^ = 2 and as{Mz) = 0.118 one obtains 
(in units of 10"^ GeV-^) 8.77 + 1.44 + 0.75 = 10.95, 8.77 -|- 1.23 + 0.72 = 10.72 and 
8.77-1-3.13-1-0.44 = 12.35 for the perturbative expansion of the moment for the first, 
second and the third method, respectively. The discrepancy between the first two 
methods and the third one is considerably larger than an error estimate based on the 
size of the order terms obtained from each of the methods. These discrepancies 
represent theoretical uncertainties that are not captured by conventional variations 
in the renormalization scale. However, they should be included in the estimate of the 
theoretical error. We include this type of theoretical uncertainty in our procedure 
to determine the combined results shown in Tab. 01 

As our final result we adopt 

rric{mc) = 1.29 ±0.07 GeV (8) 

obtained from our fits using the moments for n = 2,3. They are sufficiently below the 
duality bound and, at the same time, have only small sensitivity to the continuum 
region above the ip' resonance. Our result is compatible with earlier analyses based 
on moments of Rcc IIH ESI EH EH Ell , and with recent lattice determinations jlOl 
El 021 as well as with the finite energy sum rule analysis of Ref. jlHl- However, our 
error is much larger than in Refs. |21 EH EH] due to our more pertinent treatment of 
theoretical and experimental uncertainties. On the other hand, the size of our error 
is comparable to Refs. [SEl EHl where summations of terms proportional to powers 
of (asi/n) [HI El were partially included to account for non-relativistic effects in the 
low-energy region s ~ 4m^. However, the resulting central values are systematically 
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lower than ours by up to 100 MeV. This effect apparently arises from the summation 
of the higher order terms just mentioned and from the use of moments for n > 4. 
Both issues are problematic from the conceptual point of view. The size of our error 
is also comparable to the analysis of Ref. [IH], but our central value is about 100 MeV 
lower. This discrepancy seems to indicate that there is a systematic difference 
between the sum rules based on Eq. (0) and the finite energy moments employed in 
Ref. We note, however, that the strong influence of systematic errors in the 
derivation of the non-charm R ratio R^c and possible ambiguities in the choice of 
the upper energy cut-off represent sources of uncertainties for the finite energy sum 
rule approach that are not easy to quantify. 

Based on the numbers given in Tabs.lUandElit is straightforward to discuss how 
the uncertainty in rnc{fnc) might be further reduced in the future. One improvement 
could be achieved by the computation of the order corrections to the moments - 
if they indeed lead to a reduction of the discrepancies between "contour-improved" 
and "fixed-scale" predictions. At present the theoretical uncertainty coming from 
the different ways to implement the renormalization scale amounts to about 30 to 
35 MeV. The uncertainties coming from the J ftp e"'"e~-partial width, from the error 
in as and from the treatment of non-perturbative effects are each at the level of 
10 MeV. For the combination of these uncertainties a reduction of up to 15 MeV 
appears possible for an improved measurement of the e"'"e~-partial width and if 
smaller uncertainties in the determinations of the strong coupling are achieved. On 
the other hand, improvements in the data for ip' and for the cc continuum will not 
improve the situation significantly. 

We have also applied our method to estimate the theoretical uncertainties for 
determinations of mb(mfe) from low-ra moments of Rbf,. Using the approach described 
above we find rribi^b) = 4.22 ± 0.11 GeV. The error is only 20 MeV larger than the 
error obtained from an earlier analysis by one of the authors which was based only 
the first method with an energy-independent renormalization scale. This shows that 
the scale choice //^ ~ — has a smaller impact in the bottom quark case due to 
its larger mass. 
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